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Models  of  Molecular  Rotations 
in  the  Solid  Phase 


Fig.  1— Equipotential  lines  for  K  =  I. 


Crystal  Field  Theory: 

*  Molecular  center-of-mass 
(C.M.)  remains  fixed  at 
trapping  site  center. 

*  Rigid,  undistorted  trapping 
site. 

A.F.  Devonshire,  Proc.  Roy. 
Soc.  (London)  A153,  601 
(1936). 


Z 


X 


Rotation-Translation  Coupling: 

*  Molecular  "center-of-interaction" 
(C.l.)  occupies  trapping  site  center 
at  equilibrium. 

*  Molecular  C.M.  moves  relative  to 
trapping  site  center  during  rotation. 

M.T.  Bowers,  G.l.  Kerley,  and 
W.H.  Flygare,  J.  Chem.  Phys. 
45,3399(1966). 


rotating  molecule  trapped  in  pseudo -rotating  cage 


Pseudorotating  Cage  Model: 

*  Molecular  rotation  is  coupled  to 
synchronous  pseudorotation  of 
matrix  cage  atoms. 

J.  Manz,  J.  Am.  Chem.  Soc.  102, 
1801  (1980). 
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ortho-  and  para-hydrogen 


3  i 

Rotational  levels  of  ! 
free  Hydrogen  ! 

molecule  I 


1  12IM) 

c _ (11) 


19) 


4 


Spherical  nature 
preserved  in  solid 

Yoo 


509.9  K 


Para  Hj 
1=0 


1Jp.SK 


Metastable  impurity 
~  100  ppm  level 


Ortho  H2 
1  =  1 


I.F.  Silvera,  Rev.  Mod.  Phys.  52,  393  (1980). 
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Rapid  vapor  deposition  of  solid  pH 


M.E.  Fajardo  and  S.  Tam,  J.  Chem.  Phys.  108,  4237  (1998). 
S.  Tam  and  M.E.  Fajardo,  Rev.  Sci.  Instrum.  70,  1926  (1999). 


*  Mixed  fcc/hcp 
microstructure  in  as- 
deposited  samples;  can 
anneal  to  nearly  pure  hep. 

*  High-symmetry  trapping 
sites  in  weakly  interacting 
host  solid; 

=^>  1st  solid  matrix  host  in 
which  non-hydride  dopants 
(CO,  N2,  NO...)  are  observed 
to  rotate! 

*  Small  post-annealing 
residual  inhomogeneities; 

=^>  often  observe  sharp  IR 
absorption  lines  (-0.01  cm1); 
widths  limited  by  rotational 
dephasing/relaxation. 


61st  International  Symposium  on  Molecular  Spectroscopy,  Talk  RI03,  22  June  2006,  Ohio  State  University,  Columbus,  OH 
Approved  for  Public  Release;  Distribution  Unlimited.  AAC/PA  06-02-06-272. 


4 


(■sqe)  (|/°|)0u6o| 


12C160/pH2  IR  absorptions 


13  ppm  12C160/pH2,  d=2.9  mm. 

Assignments  from  Crystal  Field 
Theory  analysis  (T.  Momose)  with 
adjustable  coe,  Be,  De,  ae,  s2,  e3. 

"1,  ||"  =^>  polarization  vs.  hep  c-axis 
"l"  =^>  CO  in  metastable  fee  sites 

(c)  T=2.4K,  annealed. 

(b)  T=4.8K. 

(a)  T=2.4K,  as  deposited. 
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(■sqe)  (|/°|)(H6o| 


CO  isotopomers  in  pH2 


relative  wavenumber  (cm1) 


Annealed  samples  at  T  =  2.4  K. 

Spectra  shifted  by  subtracting 
gas-phase  R(0)  line. 


isotopomer 

BpH2 

(crrr1) 

Bgas 

(crrr1) 

BpH2^Bgas 

12C16Q 

1.577 

1.9313 

0.8165 

13C160 

1.538 

1.8462 

0.8331 

12C18Q 

1.469 

1.8391 

0.7988 

13C1SQ 

1.437 

1 .7540 

0.8193 

Matrix  effect  on  rotational  constant  B 
depends  on  CO  isotopic  composition. 
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pH2-gas  phase  vibrational  origin  shift  (cm'1) 


CO/pH2  vibrational  origin  shift 


gas  phase  vibrational  origin  (cm'1) 


isotopomer 

v0(pH2) 

(cnr1) 

v0(gas) 

(cm-1) 

Av0/v0 

12C160 

2140.280 

2143.254 

-0.001388 

12C170 

2113.345 

2116.279 

-0.001386 

i3Cieo 

2093.142 

2096.051 

-0.001388 

12C180 

2089.204 

2092.106 

-0.001387 

13C170 

2065.567 

2068.437 

-0.001388 

13C180 

2040.843 

2043.677 

-0.001387 

Normalized  matrix  shift: 

Av0/v0  = -0.001387(1) 
does  riot  depend  on  CO 
isotopic  composition. 
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[R(0)-P(1  )]pH2  /  [R(0)-P(1  )]s 


Beff(PH2)  I  Bavg(gas)  VS.  C.M 


0.82  I . 


0.62  0.63  0.64  0.65  0.66  0.67  0.68  0.69 

position  of  vibrationally  averaged  C.M.  (A) 

12C  C.M.  160 


0  0.65  1.14  A 


gas  phase 
isotopomer 

^avg 

(errr) 

Ravg 

(A)9 

|C-C.M.| 

(A) 

12C16Q 

1.9313 

1.1371 

0.6497 

13Q180 

1.8462 

1.1370 

0.6271 

12C18Q 

1.8391 

1.1370 

0.6822 

ISQISQ 

1.7540 

1.1368 

0.6600 

Excellent  linear  correlation  between 
matrix  effect  and  location  of  C.M. 
(plotted  relative  to  C  atom); 
r2  =  0.9963  for  fee,  0.9959  for  hep. 

But,  Rotation-Translation  Coupling 
and  Pseudorotating  Cage  models 
predict  curvature  in  Bef/B  vs.  C.M.! 
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[R(0)-P(1  )]pH2  /  [R(0)-P(1 )], 


Center  of  Interaction  for  CO/pH2 


*  Extrapolation  of  linear  correlations 
back  to  Beff/Bavg  =  1  (i.e.  vanishing  matrix 
effect  on  rotations)  yields  "experimentally 
determined"  C.l.  located  0.402  A  from  C 
atom,  or  0.248  A  from  C.M.  of  12C160. 

*  Excellent  (fortuitous?)  agreement  with 
0.24-0.27  A  C.l.  -  C.M  separation 
calculated  by  energy  minimization  in  linear 
pH2-CO-pH2  construct  for  Rnn  =  3.8-4. 1  A. 

*  Caveat!  For  N2/pH2:  C.l.  =  C.M.  and 
Bef/Bavg  =  0.90(±0.02),  riot  unity. 


S.  Moroni,  M.  Botti,  S.  DePalo,  and  A.R.W.  McKellar, 
J.  Chem.  Phys.  122,  094314  (2005). 

P.  Jankowski  and  K.  Szalwicz, 

J.  Chem.  Phys.  108,  3554  (1998). 

P.  Jankowski,  potH2CO_V04.f, 
private  communication  (2004). 
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Relevance  to  rotations  in  I  He? 


Adiabatic  Following  Regime 


1  - 


i  i  n  | - 1 - 1 — r~r~n  rq - 1 - 1 — i  rrrrq - 1 - 1 — r  i  n  n  | - 1 - 1 — i 

0.01  0.1  1  10 


Figure  11 


B  (cm"1) 

gas  '  ' 


Slightly-hindered  rotation 
regime,  where  insights  from 
studies  in  quantum  solids 
may  be  valuable. 


M.  Y.  Choi,  G.  E.  Douberly,  T.  M.  Falconer,  W.  K.  Lewis,  C.  M.  Lindsay,  J.  M.  Merritt,  P.  L.  Stiles, 
and  R.  E.  Miller,  International  Reviews  in  Physical  Chemistry  25,  15  (2006). 
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CO  /  IHe  alternative  analysis 


1.9  - 

. . 1 . . . . . 

1.80  1.85  1.90  1.95  2.00 


1  1  1  l  1  1  1 . I . 1  ■  1  1  . . |  |  ■ 

1.80  1.85  1.90  1.95  2.00 


2  H/H, 


FIG.  7.  Rotational  energies  of  CO  isotopomers  plotted  versus 
the  ratios  of  the  reduced  masses  (squares).  is  the  reduced  mass  of 
the  isotopomer  /,  with  i=  1  denoting  the  reference  12C  ,60  species. 


Assumes:  B|He/Bgas  =  constant  =  0.63 
=>  Av0  =  -0.214  to  -0.317  cm'1 

K.  von  Haeften,  S.  Rudolph,  I.  Simanovski, 

M.  Havenith,  R.E.  Zillich,  and  K.B.  Whaley, 
Phys.  Rev.  B  73,  054502  (2006). 


Assumes  Av0/v0  =  -0.0002 
(Av0  =  -0.429  cm-1  for  12C160) 

Chosen  to  give  best  straight  line 
correlation  B|He/Bgas  vs.  C.M. 
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General  applicability  ? 


Analysis  of  HCN  /  IHe 
rotational  and  ro-vibrational 
spectra  via  our  approach  is 
inconclusive. 

Need  more/better 
experimental  data  to  test 
for  linear  correlation 
between  B,He/Bgas  and 
position  of  C.M. 

^  12C  C.M  160 

| - 1- - 1- - 1 

0  1.07  1.62  2.22  A 


K.  Nauta,  Ph.D.  Dissertation,  U.  North  Carolina,  Chapel  Hill  (2000). 

A.  Conjusteau,  C.  Callegari,  I.  Reinhard,  K.K.  Lehmann,  and  G.  Scoles,  J.  Chem.  Phys.  113,  4840  (2000). 
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Quantum  Solids  vs.  Liquids 


P 


He(0 

0.0218  A-3 


1  r5  i  i  i  t  j  i  i  t  i  |  i  i  t  i  j  i  i  i  i  j~~r  i  i  ~r 


0  2  4  6  8  10 

r  (X) 


Radial  pair 
distribution  functions 
in  quantum  solids 
appear  quite 
"liquid-like"  due  to 
large  amplitude  ZPM. 


He 


(S) 


p  =  0.0335  A-3 


P^2(s) 

p  =  0.0261  A3 


r(A) 


He:  S.A.  Vitiello  and  K.E.  Schmidt,  Phys.  Rev.  B  46,  5442  (1992). 
H2:  F.  Operetto  and  F.  Pederiva,  Phys.  Rev.  B  73,  184124  (2006). 
Ne:  M.E.  Fajardo,  J.  Chem.  Phys.  88,  119  (1993). 


Ne(s) 

p  =  0.0436  A'3 


0  2  4  6 
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Summary 


*  High  resolution  IR  absorption  spectra  of  diatomic 
dopants  in  solid  pH2  provide  rigorous  tests  of  theories 
describing  molecular  rotations  in  solids. 

*  Leading  candidate  theories  cannot  explain  strong 
linear  correlation  of  Beff/B  vs.  |C.M.  -  C.l.|. 

*  Physical  significance  and  general  applicability  of 
these  correlations  still  under  investigation. 

*  Are  the  similarities  between  CO/pH2  and  CO/IHe  just 
due  to  large  amplitude  ZPM  of  nearest  neighbors? 
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